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C^ . Abstract 

o 

0> ' We present new sets of fragmentation functions for charged pions and kaons, 



both at leading and next-to-leading order. They are fitted to TPC data taken 



Q^' at energy y/s = 29 GeV and describe excellently a wealth of other e'^e'~ data on 

^S^ . charged-hadron production, ranging from ys = 5.2 GeV way up to LEP 1 energy. 



They also agree with data on the production of neutral pions and kaons, if one makes 
the natural assumption that the respective fragmentation functions are related to 
the charged counterparts by SU(2) symmetry. We also list simple parameterizations 
of the X and Q^ dependence of our results, which may be implemented conveniently 
in applications. 
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05 6 HH 93P (5) and by EEC Program Human Capital and Mobility through Network Physics at High 
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1 Introduction 

The Lagrangian of quantum chromodynamics contains quarks and gluons as elementary 
fields. Allowing for these particles to appear as asymptotic states, we can evaluate scatter- 
ing amplitudes perturbatively, in principle to arbitrary precision. Of course, this picture 
violates the principle of confinement of colour; experiments detect hadrons rather than 
quarks and gluons. Nevertheless, this simplified computational procedure is very success- 
ful in describing the production of jets of hadrons at high transverse momenta (px) in 
scattering processes and at high centre-of-mass (cm.) energies (y/s) in e+e~ annihilation. 
Apparently, clustering partons in the final state yields a valid approximation, although 
this does not account for any details of hadronization. On the other hand, experiments 
are providing us with copious information on inclusive production of single hadrons, which 
cannot be interpreted along these lines. In this case, we need a concept of describing how 
partons turn into hadrons. In the framework of the QCD-improved parton model, this is 
achieved by introducing fragmentation functions, D^{x,Q'^). Their values correspond to 



the probability that a parton, a, which is produced at short distance, Jl/Q"^, fragments 
into a hadron, h, carrying a fraction, x, of the momentum of a. Unfortunately, it is not yet 
understood how these functions can be derived from the QCD Lagrangian, in particular 
for hadrons with masses comparable to the QCD asymptotic scale parameter, A. How- 
ever, given their x dependence at some scale Qq, the evolution with Q^ may be computed 
perturbatively in QCD. This allows us to test QCD quantitatively within one experiment 
observing single hadrons at different values of pt (in the case of scattering) or taking data 
at different values of ^/s (in the case of e~^e~ annihilation). Moreover, appealing to the 
factorization theorem, we may assume that the D^{x,Q'^) functions are independent of 
the process in which they have been determined, and characterize a universal property of 
h. This enables us to make quantitative predictions for other experiments, too. 

Recently, precise data on inclusive hadron production at pp, ep, and e~^e~ colliders 
have become available. A comprehensive review of inclusive particle spectra in earlier e~^e~ 
experiments may be found in ill]. In a series of papers, two of us — partly in collaboration 



with Borzumati — have confronted pp data by CDF and UA2 and ep data by HI on 
pion and kaon production with theoretical predictions at NLO PJ^. Although we found 
reasonable agreement, we have to bear in mind that a possible limitation of our analyses 
is related to the point that the literature on pion and kaon fragmentation functions is now 
rather aged and, perhaps, slightly obsolete. The most recent sets we could locate were 
the ones by Baier et al. [§], which were updated and complemented by Anselmino et al. 
0. These sets are LO fits to low-energy e+e~ and deep-inelastic muon-nucleon data. 

Recently, NLO sets for neutral pions were presented in |^. These were fitted partly to 
e"'"e~ data in the 30 GeV range, partly to fixed-target and collider data, with 23 GeV < 
y/s < 630 GeV, and partly to data of e+e~ annihilation at 30 GeV simulated using a 
Monte Carlo event generator. Unfortunately, only the x distributions at Q^ are given in 
0; an implementation of the Q^ dependence is not available. New NLO sets have also 
been constructed for the eta meson and charged pions just recently [0, again using Monte 
Carlo data. Also here, the Q^ dependence is not given. 

Since our goal is to test quantitatively the QCD-improved parton model, we consider 
it requisite to perform an independent analysis based on real experimental information. 
Complementarily to [Q, we shall focus attention on charged pions and kaons, averaged 
separately. Similarly to 0, we analyze e+e~ data at around 30 GeV. A different av- 
enue was taken in p[, where charged-hadron fragmentation functions were extracted from 
LEP 1 data, not distinguishing between the various particle species. 

The organization of this paper is as follows. In Sect. 2, we shall introduce the formalism 
needed to extract fragmentation functions from e^e~ data on inclusive hadron production 
at NLO. In Sect. 3, we shall carry out the actual analysis and discuss our results. In Sect. 4, 
we shall check our results against data on neutral-pion and -kaon production as well as 
data on charged-hadron production without separation of the particle species. We shall 
also compare our results with those of |0]. Our conclusions will be summarized in Sect. 5. 
In the Appendix, we shall list simple parameterizations of our sets and quote their range 
of validity. 



2 Formalism 

The inclusive production of a hadron, h, by e~^e~ annihilation, 

e+e- ^{-f,Z)^h + X, (1) 

is completely characterized by three observables. These are the energy fraction of the 
outgoing hadron, 

" = 7i ' ^^ 

its angle with the beam axis, 6, and the total energy of the system, ^/s. For massless 
hadrons, definition (Q) agrees with the convention used in experiments, whereas for non- 
vanishing masses, rrih, one has to adjust the theoretical prediction by including a factor 



j3h = J I — mfjE\. Choice of reference frame is not an issue here, since, at e^e~ colliders, 
the laboratory and cm. frames coincide. The cross section of (|I]) exhibits the angular 
structure 

dxdcosO 8 dx 4 dx 4 dx ^ 

where the superscripts T and L denote the contributions due to transverse and longitu- 
dinal polarizations, respectively, and the asymmetric term, labelled A, accounts for the 
interference of the photon with the Z boson. Usually, experiments do not determine 6 
distributions. Thus, we shall integrate over 6. This eliminates the asymmetric term in 

(0)- 

The partonic subprocesses may be treated exactly in perturbative QCD. However, this 

is not yet possible for the process of hadronization. As explained in the Introduction, the 

latter is described in terms of phenomenological fragmentation functions, which must be 

extracted from experiment. In the language of the QCD-improved parton model, the x 

distribution of (|l]) emerges from the x distribution of e"'"e~ — > a + X, {daa/dx){x, /x^, Q"^), 

through convolution with D^{x^ Q^), 

s da{e~^e~ —>■ h + X) .^-^ f ^^ r^h/- i\j2\ daafx 



/3h dx a^z ^ dy 



^/^Z,;(..M|).^(£.,^M|). (4) 



Here, the sum extents over all active partons (a = g,u,d,s,c,b), fi is the renormaliza- 
tion scale of the partonic subprocess, and Mf is the so-called fragmentation scale, which 
corresponds to the typical energy where the transition from the perturbative to the non- 
perturbative regime takes place. In next-to- leading order (NLO), Mf defines the point 
where the divergence associated with coUinear radiation off parton a is to be subtracted. 
As usual, we work in the massless-quark approximation, i.e., we neglect threshold effects. 
As a consequence, the bottom quark is taken to be active throughout the entire Q^ evo- 
lution and to contribute to the cross section also at low energies. We checked that this 
slight inconsistency is inconsequential for the goodness of our fits. To NLO in the MS 
scheme, the cross sections of the relevant subprocesses are given by M 



da, 



dy 
da 



-{y,fx\MJ)=a^N, 



5{l~y) + 



asifi"^ 
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'2Nf 



iy,fi',Mf)=aoNjY.< 



2n 



nrn2/)ln- 



.. v.."^^T + ^J(2/) + ^."(y) 



(5) 



where (Tq = (47r«QEj~|/3s) is the total cross section of e"'"e~ — >• /i"*"//" for massless leptons 



(0,T), 



Nc = 3, Cq. is the electric charge of quark gj in units of the positron charge, P^^' (y) is 



the LO term of the transverse a —^ b splitting function. 



Pa,iy,c^s{Q')) = Pr'{y) + 



(0,T), 



2n 



,(i,r) 



P^,'"{y) + 



(6) 



and 



Kliv) = Cf 



2^ 



:i + y' 
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S{l-y) 



K^{y) = Cf, 



Kiy) 



c. 
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1-2/ 



\\i{l-y) + 2\ny] -2 
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(7) 



with Cp = (A'c — l)/(2A'c) = 4/3. As usual, the subscript -|- denotes distributions. The 
NLO formula of as(/i) may be found, e.g., in [|l^. Throughout this work, we choose 



A = 190 MeV, in order to be compatible with recent sets of parton density functions for 
the proton and photon. The scales fi and Mf are identified with the only intrinsic energy 
scale, ^/s. For this choice of scales, the logarithmic terms in (^) involving s vanish, so 
that the NLO corrections are expressed just in terms of the K functions. We observe that 
^Jiy) + ^g{y) < fo^ ^^y value of y, so that the gluon contributes destructively to @j. 
Having defined the partonic subprocesses, we turn to the fragmentation functions. 
Their x distributions are not yet calculable in the framework of perturbative QCD. How- 
ever, once we know them at some scale Qq, the Q^ evolution is determined by the DGLAP 



equations [TT[]. Our task is thus to construct a model for the x distributions at a starting 
scale Qo! which, after evolution, fits the data at scale Q^ . The DGLAP equations read 

That is, we have to solve a system of integro-differential equations. This may be achieved 

with the help of the Mellin-transform technique |12[. The Mellin transform, F, of a 

function, /, is defined as the set of its moments, 

1 
F{n)= I dxx''-^f{x). (9) 



The inverse transformation is 

7+joo 
7— ioo 

where the path of integration must be chosen to the right of all singularities of the in- 
tegrand. In our numerical analysis, we choose 7 = 10. The essential property of this 
transformation is that it renders convolutions to products. In fact, denoting the moments 
of P^ and D'^ by Aah and M^, respectively, we have 



^ ^4-<^ r^2^ _ <^s{Q 



2\ 






^^Mcin, Q') = ^^^ {A,,{n, as{Q'))M,,{n, Q^) + A,,(n, as{Q^))MG{n, Q^)) , (11) 



where A^s refers to the usual non-singlet combination of quark and antiquark splitting 



functions fl2 , and 



Nf 

MG{n,Q') = Mg{n,Q'). (12) 

The LO expressions for Aab may be found in [|^,|^[] and the NLO expressions in [0]|. 
Note that the non-singlet terms, Mf, decouple from the gluon. The solutions of ([TI|) 
factorize, 

Min,Q') = EiQl,Q')Min,Ql). (13) 

The relevant evolution operators, E^Q^jQ"^), may be found in |]12|. We parameterize the 
X dependence of the fragmentation functions at Ql as 

D{x,Ql) = Nx''{l-xf. (14) 

The corresponding moments read 

M{n,Ql) = NB{n + a,f3 + l), (15) 

where B is Euler's beta function. The evolution in n space is performed analytically via 
(p!3|). Finally, the D{x,Q'^) functions are obtained by solving (|10D numerically. This re- 



quires knowledge of Aab for complex n. A prescription for the proper analytic continuation 
may be found in |T3|. 

In (p^ , we introduced three parameters for each parton. To reduce the number of free 
parameters, we shall discriminate only between valence- and sea-type quarks. To that 

^Note that there is a typographical error in (B.46c) of [Q. Cp should be replaced by Ca- 

^Note that there is a typographical error in (A. 10) of ||lj|. The term before the last one of the 

contribution proportional to CfTrJ should read — (16/3(n + 2)^) instead of — (16/3(n + 1)^)- 



end, we make the following ansatz for the average of charged pions, 

Dl = D: = Dl = D: = Dl = Dl = Dl (16) 

For the average of the charged kaons, we set 

D§ = D^ = Df = Df = D^ = Df = Df . (17) 

Furthermore, we put a = —1 for valence and sea quarks, which may be motivated by 
inspection of the energy dependence of hadronic multiplicities. In the case of the gluon, 
we found that the choice a = is favoured by the data. We are thus left with six 
parameters to be determined from experiment. 

3 Results 

In general, data of inclusive hadron production at e^e^ colliders are most suitable for the 
extraction of fragmentation functions; in the case of fixed-target, collider, and ep data, 
the information on the fragmentation functions is obscured by theoretical uncertainties 
arising from the parton density functions and the choice of factorization scales connected 
with the initial state. For our analysis, we select the data on charged-pion and -kaon 



production taken at energy ^/s = 29 GeV by the TPC Collaboration at SLAC |jT6|. These 
data combine small statistical errors with fine binning in x and are thus more constraining 
than data collected by other e~^e~ experiments in the energy range from 5.2 to 44 GeV. 
For this reason, our approach is to fit exclusively to TPC data and to use the other data 
for cross checks. 

The quality of the fit is measured in terms of the average Xd o f ^^ ^11 data points 
with X > Xcut = 2 GeV/y/s. This restriction is necessary to exclude events in the non- 
perturbative region, where the above formalism is bound to fail. Our technical procedure is 
as follows. We consider (§) at /i = Mj = y/s = 29 GeV as a function of the six parameters 
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in LO and 



that determine the x dependence of the fragmentation functions at Q^ = 2 GeV^. Using a 
multidimensional minimization algorithm [|17|, we search this six-dimensional parameter 
space for the point at which the deviation of the theoretical prediction from the data 
becomes minimal. Our results are listed below. For the average of charged pions, we 
obtain 

D^y'''^\x, Qo) = 0.551 x-^ (1 - x)^•^ 
dJ '"-""^(x, Qo) = 1.23 a:-^ (1 - x)^-'\ 
4^'^°)(x,Qo)=3.77(l-x)2-° (18) 

jji-^NLO)^^^ Qo) = 0.338 x-^ (1 - xy-^ 

4^'^^^)(x, Qo) = 1.19^-^1 - a:)'■'^ 

Dj:'^^°)(a;,Qo) = 3.65(l-a;)2-o (19) 

in NLO. In the case of kaons, we find 

Djf '^°^(x, Qo) = 0.155X-1 (1 - xy-\ 

DP°)(a;, Qo) = 0.216a:-i (1 - x)l■6^ 

Dg"'^°^(x, Qo) = 0.324 (1 - x)2-° (20) 

in LO and 



<'""^^(x,Qo)^ 


= 0.099x-^(l-x)^■^ 


4"'"'^^(-,Qo)^ 


= 0.214x"^(l-x)2■^^ 


4"'"'^^(a^,Qo)^ 


= 0.372 (l-x)2-° 



(21) 

in NLO. For the reader's convenience, we list simple parameterizations of the x and 
Q^ dependences of these sets in the Appendix. It cannot be expected that such simple 
parameterizations reproduce the true result perfectly for all x and Q^ values that occur 
in applications. Deviations in excess of 10% may occur at Q^ < 25 GeV^ and x < 0.1, 

9 



in particular for the gluon. Nevertheless, we believe that such parameterizations are 
indispensable for practical purposes, especially at NLO. 

Our fragmentation functions exhibit good perturbative stability. In the case of the 
pions, there is only little difference between our LO and NLO sets. The difference is 
somewhat bigger for the kaons. In the comparisons with data, the NLO sets do just 
slightly better. The x dependence of (|T9| , ^TD at Ql is visualized in Figs. |^J|, respectively. 
In the following, we shall always employ the parameterizations given in the Appendix. 
The goodness of our fits to the TPC data on charged-pion and -kaon production may 
be judged from Figs. 0,|[ respectively. In LO (NLO), we achieve xio.f. — 0-'^ i^-'^) ^^ 
pions and 0.6 (0.5) for kaons. Here and in the following, Xdof refers to the interval 
max(xcut; 0.1) < X < 0.9. Notice that our parameterizations fit the TPC data excellently 
even down to x = 0.06. 

In Figs. |-0, we postdict data taken by DASP ^ and ARGUS [^ at y/s = 5.2 and 



10 GeV, respectively. We have made comparisons also with TASSO data at ^/s = 34 and 
44 GeV pO|! which we do not, however, display graphically. Data of inclusive charged- 



kaon production are well reproduced by our fits. In Fig. |^, this is nicely demonstrated for 
ARGUS. In the case of the charged-kaon data by DASP, for which we do not show a figure, 
we reach Xdot — 1-4 (1-2) in LO (NLO), which we can bring down to values below unity 
when we exploit the freedom in the overall normalization of the data by about ±12%. 
The charged-kaon data by the other experiments |T0,[TP|,|2D| yield Xdof values between 0.4 



and 0.6 without adjustment, both in LO and NLO. 

The situation is slightly less favourable in the case of charged pions. The pion data by 

DASP [l^ behave similarly as their kaon data (see Fig. ^). However, the ARGUS data 

T9| are inconsistent intrinsically at large x, which leads to increased Xdof values, 2.2 in 



LO and 1.8 in NLO (see Fig. ^). For x < 0.1, the theoretical curves overshoot the data. 
This is presumably attributed to mass effects, which have been ignored in our approach. 
As per construction, the TPC data are reproduced well, with Xdof — ^-"^ both in LO 



and NLO. The TASSO data ||2^ are fairly described, with Xdof values between 1.4 and 



2.0. However, these data are of limited use for our purposes because, for x > 0.2, their 
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binning is a factor of 3 to 5 coarser than that of the TPC data. Furthermore, the data 
sets at y/s = 34 and 44 GeV do not comply with the expectation that scahng violation 
renders spectra at higher energies steeper. Leaving aside the TASSO data, the goodness 
of our charged-pion fits is quite satisfactory. 

4 Applications 

Since fragmentation functions are assumed to be universal, i.e., independent of the process 
from which they have been extracted, we can apply them to make predictions for inclusive 
particle production at fixed-target, collider, and ep scattering experiments, too. In the 
following, we shall confine ourselves to more extended studies of e~^e~ annihilation, which 
will substantiate our confidence in the virtue of our sets. Work on inclusive particle 
production in NLO at HERA and pp colliders is in progress. 

We commence by elaborating a welcome generalization of our results. By SU(2) sym- 
metry, our charged-pion sets are valid also for the neutrals. This observation is supported 
by the inspection of neutral-pion data, e.g., by ARGUS E^, TPC P[, JADE pi, and 



CELLO p4[, which agree with the respective charged-pion data within the errors quoted. 
This is reflected also by the Xd o f analysis. Cross checks with new data by L3 I^S] prove 
to be very successful, too. By analogous theoretical considerations, we can establish 
a relationship between the fragmentation functions for the average of the K^ and K 
mesons and those for the charged kaons. Specifically, the u, u and d, d fragmentation 
functions swop places in (p!7|). Again, we find reasonable agreement way up to LEP 1 



energy [|I|,0,0,g3J2l 

As a second application, we deal with inclusive charged-hadron production. I.e., we 
consider experiments that observe charged tracks without identifying the particle species. 
In this context, it has become customary to sum over the produced hadrons rather than 
to form the average. Clearly, pions and kaons make up the bulk of the cross section. 
By including also protons and antiprotons, the next important subset, we obtain a good 
approximation for the sum of all charged hadrons. Inspection of TPC data [|l^ reveals 



11 



that the averaged cross section of inclusive production of protons or antiprotons has a 
shape similar to the one of charged pions, but is reduced in size by a factor of 0.16 ± 0.02. 
Thus, we take the hadronic cross section to be 



dx 



1 + 0.16) + 



(22) 



dx dx 

We compare this prediction with high-precision data by DELPHI [^ in Fig. || and find 
good agreement, with xl.o.f. — O-^ (1-0) i^ LO (NLO). Our results also agree well with data 



by MARK II [g, AMY [0, TASSO fS^, and CELLO Q. Whereas the x dependence 
is reproduced perfectly for the MARK II, AMY, and DELPHI data, there are minor 
deviations in the case of TASSO and CELLO. Except for DELPHI, the Xdof values 
exceed unity in all these cases, which indicates that the incorporation of baryons by an 
overall Q^-independent factor is perhaps to naive. 

In order to test scaling violation, we compute (p2D at y/s = 91 and 29 GeV, plot the 
ratio as a function of x in Fig. 0, and compare the answer with the corresponding ratio 



of DELPHI [29| and MARK II data ^^, taking into account only the statistical errors 



of the latter. For x > 0.1, our prediction agrees remarkably well with the experimental 
result, both in LO (dashed line) and NLO (solid line). We note in passing that a similar 
analysis was carried out recently by Nason and Webber |^, who encountered a certain 
mismatch between their calculation and their choice of data. These authors determined a 
set of fragmentation functions by fitting to the DELPHI data on inclusive charged-hadron 
production |^9|, evolved their results backwards down to PETRA energy, and compared 



the outcome with TASSO data [^, taking into account only the statistical errors of 
the latter. They interpreted the discrepancy that they discovered by suggesting that 
significant nonperturbative power corrections are necessary to reproduce the measured 
scaling violation (see Fig. 14 of P|). In our opinion, the major source of this problem is 
the neglect of the systematic errors on the TASSO data rather than the influence of power 
corrections. In fact, when we repeat the analysis of Fig. ^ for TASSO and CELLO [Q 



data, both taken at PETRA, again including only the statistical errors in the denominator, 
we find that the two sets do not overlap and are separated by our curves. This suggests 

12 



that the TASSO and CELLO data are not consistent with each other when only their 
statistical errors are taken into account. 

Another look at scaling violation is presented in Fig. |10|, which shows the Q^ depen- 



dence of the charged-hadron cross section for various values of x. The evaluations of (22) 



to LO (dashed lines) and NLO (solid lines) are confronted with a rich collection of data 



lq , p!9| , |29| , |30| , |3Ti|32| , p3| , which have been interpolated in x where appropriate. We conclude 
that the data indeed exhibit the expected scaling violation. 

The DELPHI Collaboration used an alternative approach to interpret their data |2^ . 



The inclusive hadron spectra are generated from a Monte Carlo program based on order 
a| QCD matrix elements, which produces jet cross sections and takes into account the 
subsequent fragmentation of the jets into hadrons according to certain model assumptions. 
In this approach, the cancellation of infrared and coUinear singularities occurs through the 
definition of jet cross sections with chosen jet resolution criteria. This deviates from the 
canonical factorization approach, where the coUinear singularities are absorbed into the 
fragmentation functions. We strongly advocate the canonical formalism, which is intrinsic 
to the QCD-improved parton model. 

At this point, we should address a possible shortcoming of our analysis, which is 
common to all studies of fragmentation in e^e^ annihilation, namely the circumstance 
that the gluon fragmentation function is poorly determined. This is due to the fact that 
the gluon participates in the process only at NLO, and even then has an appreciable 
impact on the cross section only at very small x values close to Xcut- Thus, it contributes 
mainly through the Q^ evolution, where it couples to the singlet combination of quarks. 
Consequently, the gluon parameters are strongly correlated with the sea-quark parameters 
and are feebly constrained by the data. Hadron and ep scattering experiments deliver 
far more information on the gluon fragmentation function. Fortunately, the new three-jet 
data by OPAL ||3^ are able to probe gluon fragmentation directly. In Fig. |TT], we compare 



these data with our LO (dashed line) and NLO predictions (solid line) according to (|22D . 
The agreement is reasonable. We tried also different ansatze for Dq{x^Q'^^ by choosing 
a 7^ in (ITW). While they worked similarly well for all e^e~ data considered above, they 
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TABLE I. Left-hand side of (|3|) at NLO for VQ^ = 5, 29, and 91 GeV. We use 
min(xcut, 0.1) as the lower bound of integration and sum over pions, kaons, protons, and 
antiprotons. 



a 


VQ^ [GeV] 


5 


29 


91 


u, d 

s 

c,b 

9 


0.81 
0.55 
0.45 
0.62 


0.76 
0.52 
0.45 
0.58 


0.94 
0.68 
0.60 
0.83 



failed in the case of Fig. |TT| p5[] . 

Finally, we investigate if our fragmentation functions satisfy the momentum sum rules. 

Guided by the idea that a given outgoing parton, a, will fragment with 100% likelihood 

into some hadron, /i, and that momentum is conserved during the fragmentation process, 

we expect that 

Y^f'dxxD^^{x,Q') = l (23) 

h J^ 

holds for any value of Q^. In our analysis, the left-hand-side of (^) should be smaller 
than unity, since we consider only pions, kaons, protons, and antiprotons, while we do not 
include hyperon channels. Furthermore, we are forced to introduce a finite lower bound 



of integration in (|23|), since the x — > limit of D^{x,Q ) is beyond our control. a;min is 
a plausible choice for this cutoff, except at small Q^ {VO^ < 20 GeV), where we use 0.1 
instead. In Table I, we list the results obtained with our NLO fragmentation functions 
for \/Q^ = 5, 29, and 91 GeV. In fact, all entries are smaller than one as they should. 



5 Conclusions and Outlook 

We presented LO and NLO fragmentation functions for charged pions and kaons in the MS 
scheme. They were extracted from e"'"e~ data on inclusive production of these particles 
collected at ^/s = 29 GeV. We checked that these fragmentation functions are in agree- 
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ment with the majority of the e"'"e~ data in the energy range between 5.2 GeV and Mz-, 
and suggested possible explanations, where this was not the case. This work is motivated 
by the circumstance that a great amount of high-quality data on inclusive hadron produc- 
tion has been accumulated in fixed-target, hadron collider, ep scattering, and e^e~ anni- 
hilation experiments since the early 1980's, which was not accompanied by any progress 
in the phenomenology of charged-particle fragmentation functions. The quest for such 
progress has been increased dramatically by the advent of new high-statistics data from 
HERA and LEP 1. In particular, meaningful tests of the QCD-improved parton model 
can be performed only beyond LO, so that knowledge of fragmentation functions with 
NLO evolution is essential. The present paper makes an effort to fill this gap. 

By the factorization theorem, the use of our fragmentation functions is not restricted 
to e'^e~ collisions. They may be applied to any other process by which charged pions 
and kaons are produced inclusively. However, one should bear in mind that our sets 
come with a certain range of validity, which is specified in the Appendix. For instance, 
the production of low-p-p hadrons in high-energy hadron-hadron and ep scattering has a 
tendency to probe small values of x that lie outside this range. Moreover, such hadrons 
are likely to originate from gluons, whose fragmentation function is not yet constrained 
so well by e^e^ data. This leaves room for further improvements in the future. 
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A Parameterizations 

For the reader's convenience, we present here simple parameterizations of the x and Q^ 
dependence of our fragmentation functions. It is necessary to quote a range of validity for 
these parameterizations. In the case of Q^, they agree with the results obtained by explicit 
evolution to better than 10% if 5 GeV < v^Q^ < 91 GeV. As for x, we are experiencing 
technical limitations connected with the stability of the numerical integration for a; — > 
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and X — i> 1. Our results should be reliable in the interval 0.1 < x < 0.9. The upper bound 
on X may be ignored in practice, since the fragmentation functions take very small values 
at X = 0.9 and approach zero for x -^ 1. 

As usual, we introduce the scaling variable 

ln(Q§/A=) ' ' 

and parameterize our results by simple functions in x with coefficients that are polynomials 
in s. We find that the template 

D{x,Q^)=Nx''{l-xf{l + xy (25) 

works best in our case. We now list the parameters to be inserted in (|25|) for charged 
pions and kaons both in LO and NLO. At Q^ = Qq, these parameterizations coincide with 



Iq-pT]) as they should. 
1. LO fragmentation functions for (7r+ + 7r^)/2: 



AT = 0.551- 0.053s -0.032s2 
a = -1.0 -0.0272s 
/3 = 1.2 + 0.67s 
7 = -0.393s (26) 



4-"'^)(x,Q2) 



A^ = 1.23 + 2.85s -1.60s^ 



a = -1.0 + 0.447s -0.266s2 



/3 = 4.77- 2.88s + 2.055^ 

7 = -9.01s + 4.36s2 (27) 
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4"^^)(a:,g2): 



AT = 3.77- 5.49s + 2. 14s^ 



a = -1.73s + 0.329s^ 



f3 = 2.0 + 1.00s + 0A4f 

-f = -2.05s + 2ASf (2^ 



2. NLO fragmentation functions for (tt^ + vr )/2: 



Di;'^^^)(x,Q2): 



4^'^^^)(x,g^ 



4^'^^^Hx,Q2) 



iV = 0.338 - 0.064s - O.OlOSs^ 
a = -1.0 -0.059s 
/3= 1.2 + 0.60s 
7 = -0.163s (29) 



A^ = 1.19 + 4.20s -2.865^ 



a = -1.0 + 0.757s -0.537s2 
/3 = 5.26- 5.22s + 3.62s2 

7 7,2 



7 = -13.6s + 8.17r (30) 



AT = 3.65- 6.00s + 2.52s 



a = -1.49s-0.213s2 



-2 



/5 = 2.0- 0.23s +1.525^ 

7 = -6.53s + 8.43s2 (31) 



3. LO fragmentation functions for {K"^ + K )/2: 
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Df''^\x,Q^) 



AT = 0.155- 0.010s -O.OlOs^ 
a = -1.0 -0.012s 
/5 = 1.4 + 0.66s 
7 = -0.572s (32) 



Df''^\x,Q-) 



4"''^Ha:,Q^) 



A^ = 0.216 - 0.0744s 
a = -1.0 -0.196s 
/5 = 1.69 + 0.69s 
7 = _0.195s (33) 

N = 0.324 - 0.506s + 0.2058^ 
a = -1.860s 
/3 = 2.0 + 1.67s 
7 = -0.235s + 2.68s2 (34) 



4. NLO fragmentation functions for {K~^ + K )/2: 



Df'^^^)(x,Q2): 



AT = 0.099 - 0.0198s - 0.00275^ 
a = -1.0 -0.060s 
/9= 1.4 + 0.60s 
7 = -0.233s (35) 

N = 0.214 - 0.00073s - 0.00635^ 
18 



a = -1.0 -0.0067s 

/9 = 2.16 + 0.57s 

7 = -0.912s (36) 



Djf'^^^)(x,Q2): 



N = 0.372 - 0.606s + 0.2535^ 



a = -2.55s + 0.584s2 



/5 = 2.0- 0.52s + 1.82s2 

7 = -7.43s + 9.24s^ (37) 
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FIGURE CAPTIONS 

Figure 1: x distributions of the NLO pion fragmentation functions at Qq = 2 GeV^ 
Figure 2: x distributions of the NLO kaon fragmentation functions at Qq = 2 GeV^ 
Figure 3: Cross section of inclusive charged-pion production at ^/s = 29 GeV as a function 



of X. The LO (dashed hne) and NLO fits (sohd hne) are compared with TPC data |T6 



Figure 4: Cross section of inclusive charged-kaon production at -^i = 29 GeV as a function 



of X. The LO (dashed line) and NLO fits (solid line) are compared with TPC data |[16 
Figure 5: Same as in Fig. | for v^ = 5.2 GeV and DASP data [0] 



Figure 6: Same as in Fig. for y^ = 10 GeV and ARGUS data ||T9 



Figure 7: Same as in Fig. | for v^ = 10 GeV and ARGUS data p 

Figure 8: Cross section of inclusive charged-hadron production at ^/s = 91 GeV as a 
function of x. The dashed (solid) line corresponds to the prediction on the basis of (|22|) 
at LO (NLO) 

Figure 9: Ratio of the x distributions of inclusive charged-hadron production at ^/s = 91 



and 29 GeV. The dashed (solid) line corresponds to the prediction on the basis of (|2^ at 
LO (NLO). The points correspond to DELPHI data [g|] divided by MARK II data fS^, 
taking into account only the statistical errors of the latter 

Figure 10: Scaling violation in the cross section of inclusive charged-hadron production. 
Shown is the Q^ dependence for x = 0.075,0.15,0.25,0.35,0.5,0.75,0.9. The dashed 
(solid) line corresponds to the prediction on the basis of (^) to LO (NLO). The data are 
taken from [|18|,|19],^,^,0,^,Q and are partly interpolated in x 

Figure 11: Cross section of inclusive charged-hadron production at ^i = 91 GeV as a 
function of x. The dashed (solid) line corresponds to the prediction on the basis of (^) at 
LO (NLO). The data are taken by OPAL |^ and represent the first direct measurement 
of the gluon fragmentation function for charged hadrons 
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